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Abstract 

Supersymmetric loop contributions can lead to different decay rates of H + — > tb and 
H~ — > bt. We calculate the decay rate asymmetry 8 = [T(H + — > tb) — T(H~ — > 
bt)] I [T(H + -> tb) + T(H~ -> bt)] at next-to-leading order in the MSSM with complex 
parameters. We analyse the parameter dependence of S CP with emphasis on the phases 
of At and A^. It turns out that the most important contribution comes from the loop 
with stop, sbottom, and gluino. If this contribution is present, 8 can go up to 10— 15% 
for tan (3 ~ 10, and to ~ 5% for large values of tan (3. 
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1 Introduction 



Already for some time, it has been customary to look for CP-violating effects beyond the Stan- 
dard Model (SM). All extensions of the SM contain possible new sources of CP violation through 
additional CP-violating phases. In particular, in the Minimal Supersymmetric Standard Model 
(MSSM), the Higgs mixing parameter fi in the superpotential, two of the soft SUSY-breaking 
Majorana gaugino masses Mj (i = 1,2,3), and the trilinear couplings Aj (corresponding to 
a fermion /) can have physical phases, which cannot be rotated away without introducing 
phases in other couplings Q. From the point of view of baryogenesis, one might hope that 
these phases are large 0] . On the other hand, the experimental limits on electron and neutron 
electric dipole moments (EDMs) §, \d e \ < 2.15 x 1(T 13 e/GeV, \d n \ < 5.5 x 1(T 12 e/GeV, 
place severe constraints on the phase of /i, <ft^ < O(10~ 2 ) M, for a typical SUSY mass scale 
of the order of a few hundred GeV. A larger M imposes fine-tuned relationships between this 
phase and other SUSY parameters ||. CP- violating effects that might arise from A u ^ (where 
u, d are light quarks) are very much suppressed as they are proportional to m u> d- On the other 
hand, the trilinear couplings of the third generation A t fi, T can lead to significant CP- violation 
effects, especially in top quark physics || . Phases of \i and A t ^ T also affect the Higgs sector in 
a relevant way. Although the Higgs potential of the MSSM is invariant under CP at tree level, 
at loop level CP is sizeably violated by complex couplings |7], §. As a consequence, the three 
neutral Higgs mass eigenstates are superpositions of the CP eigenstates h°, H°, and A . 

In this paper, we study CP violation in the decays H + — > tb and H~~ — > bi in the MSSM 
with complex parameters. In particular, we calculate the CP-violating asymmetry 

CP _ r (g+ -> tb) — r (h~ — » bt) 

r (h+ tb) + r (h- -> bt) ' u 

due to one-loop exchanges oft, b, g, x^, X°; an d H°, see Fig. 1, taking into account CP violation 
in the neutral Higgs system according to ||. Of course, the diagrams of Fig. 1 only contribute 
to 8 if they have an absorptive part. Since M is highly constrained, the most important 
phases in our analysis are <f> t and fa, the phases of A t and A\>. Therefore, we expect the graph 
with t, b, and g in the loop to be the most important one, and 8 CP to be large in the case 
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m#+ > mfj + m bi . In principle, there would also be a contribution due to v and f exchange 
analogous to Fig. le. However, this can be neglected in our study. 

The paper is organized as follows: In Section 2 we give the basic formulae for the H ± — > tb 
decays and define the decay rate asymmetry 8 CP at the 1-loop level in terms of CP-violationg 
form factors 8Yf p (i = t,b). The explicit expressions for SY t c b p due to the diagrams of Fig. 1 
are given in Section 3. In Section 4, we perform a detailed numerical analysis. In Section 5, we 
summarize our results and comment on the feasibility of measuring the CP- violating asymmetry 
8 . Appendices A, B, and C contain the necessary mass and mixing matrices, the couplings, 
and the definition of the two- and three-point functions used in this paper. 

2 The if* decay 

The matrix elements of the H + — > tb and H — > bi decays can be written as 

M H+ = u(p t )[Y+P R + Y+P L ]v(-p b ), (2) 
M H - = u(p b )[Y t -p R + Y b -P L ]v(-p t ), (3) 

with P RL = |(1 ±75) and the loop-corrected couplings 

Y± = y i + 8Y+ (i = t,b)\ (4) 

y t and yb are the tree-level couplings, 

y t = h t cos £}, yb = h b sin (3, (5) 

with h t and h b the top and bottom Yukawa couplings. The decay widths at tree level are given 
by 

r ° (Jf± _> tb) = 16?r 3 ^3 [(m 2 H+ - m\ - m 2 b ){y 2 t + y 2 b ) - Am t m b y t y b ] , (6) 

where k = K(m 2 H+ ,m 2 , ml), k(x, y, z) = [(x — y — z) 2 — Ayz] 1 ^ 2 . Since there is no CP violation 
at tree level, T°(H + -> tb) = T°(H~ -* bi). At next-to-leading order (NLO) we have 

r (iy± _ tb) = 167r 3 ^ 3 [{m 2 H+ - m 2 t - m 2 b )(y 2 + y 2 b + 2y t Re 8Y± + 2y b Re 8Y±) 

- Am t m b {y t y b + y t Re 5Y± + y b Re 6Y±j\ . (7) 
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The form factors SYf 1 (i = t, b) have, in general, both CP-invariant and CP-violating contribu- 
tions: 

5Y± = 5Yl nv ± \ SY^ . (8) 

Both the CP-invariant and the CP-violating contributions have real and imaginary (absorptive) 
parts. CP invariance implies that the form factors of H + and H~ are equal. Using Eqs. (|7]) 
and (H), we can write the CP-violating asymmetry S CP of Eq. ([!]) as 

5 cp = A (y t Re 6Y t CP + y b Re 6Y b CP ) - 2m t m b (y t Re 6Y b CP + y b Re 5Y t CP ) 

&{y? + vl) -4m t m b y t y b 

where A = m 2 H+ — mf — m b . 

At one loop, there are six generic diagrams that may contribute to 5 CP . These are (A) tri- 
angle diagrams with (i) two fermions and a scalar, (ii) two scalars and a fermion, and (hi) a 
scalar, a vector boson and a fermion in the loop, and (B) H + -W + self-energy diagrams with 
(i) two fermions, (ii) two scalars, and (iii) a scalar and a vector boson in the loop. In the 
following, we work out the formulae for 5Y^ P for the specific case of the MSSM with complex 
phases. The relevant Feynman diagrams are shown in Fig. [TJ. 



3 CP-violating contributions 

3.1 Generic diagrams 

According to Eqs. (Q) and (^ we write the matrix elements of the 1-loop diagrams of Fig. [T] as 

M H+ = u( Pt ) [5Y b + P R + 6Y+P L ] v(- Pb ) , (10) 

and analogously for the H~ decay. In fact, we are only interested in the CP- violating parts 
Re 5Yf p . Since Re 5Yf p = Re (5Y t + — 5Y~), we just need the structure lm(g gig 2 ) x Im(PaVe) 
of the form factors 5Y b + and SY t + . Here gog\g2 stands for the product of the couplings and PaVe 
for the Passarino-Veltman two- and three-point functions B Q and Cq^^- in the following, 
we give the formulae for the various contributions to Re5Y t c b p . The necessary MSSM mass and 
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Figure 1: Sources for CP violation in H + — > tb decays at 1-loop level in the MSSM with complex 
couplings (i,j = 1, 2; k = 1, 4; / = 1, 2, 3). 

mixing matrices, the couplings, as well as the definition of the two- and three-point functions 
are given in Appendices A, B, and C. 



3.2 Vertex graphs 

3.2.1 Neutralino— chargino stop (sbottom) loop 

The graph of Fig. la, with a neutralino, a chargino, and a stop in the loop, leads to 
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m t m b Im(i^yS) + m x +m x° k + ™t™*t M*$*&*£ 



+ m b m tl Im(F^ayg) + m| lm(F^ k a\ k k 



i J ' 



+ m t 
+ m h 



m Jm(F/ fe 4 fc A;g) + m^o Mif^J) + m 6 Im(i^&yg) + Im^X^ 



mt Im(F^4Zg) + m^o Im(F^a* fc /g) + m 6 Im(F/ fc a* fe fc**) + m~+ Im(F/ fc a* fc /g; 



+ M^<4*&) Im(5 (m^ +) m|o, m|+)) 



Im(Co) 
Im^) 
Im(C 2 ) 
(11) 



with Cx = Cxijnl, m 2 H+ , m\, m~ , m| , m ~+); -X" = 0,1,2, the three-point functions in the 



notation of |TTJ; Re 8Y t CP (x k xfti) is obtained from Eq. ([□]) by interchanging 
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(Ku <r+ b\ 



The contribution from the neutralino-chargino-sbottom loop has exactly the same structure. 
Therefore, Re SY b CP (x^Xj h) 1S obtained from Eq. (|ll]) by the following substitutions: for the 
masses of the loop particles m^o — > rn^+, m-+ — > m.^, rrif. — > m^; for the couplings a\ k — > l^, 
b\ k -> fc*., /eg -> &£, and Zg -> a&; and analogously for Re5F 4 CP (x°X^i)- 

3.2.2 Stop— sbottom— neutralino loop 

The stop-sbottom-neutralino loop of Fig. lb gives 

ReSY^ikbjxl) = ^ [m^m(G 4tJ aib%)lm(C ) 

-m t \m{G^b\ k b%) Im(d) 
He^ftft^) = ^ [m 4 Im(G 4 ^4a5* fc )Im(C ) 

— m 4 Im(G4^a- fc a^) Im(Ci) 

with C x = C x (m? , m 2 H+ , , m|o , m? , ra| ) . 

3.2.3 Stop sbottom— gluino loop 

The contribution from the diagram with a stop, a sbottom, and a gluino in Fig. lb is 

Re5Y b CP (U~bfg) = [m- g Im^A*^ e*>) Im(C ) 

+ m t ^(G^i?!*^) Im(Ci) + m b Im (G^i^i^) Im (C 2 ) ] > ( 14 ) 

~ 4 a - 

O 7T 

+ m* Im^X*^) Md) + m b Im^i^f,) Im(C 2 ) ] , (15) 
with C* x = C x (m 2 t ,m 2 H+ ,m 2 b ,m 2 -,m 2 k ,ml ) and a s = g*/(4m). 



- m fc Im(G 4jJ aL< fc ) Im(C 2 ) ] , (12) 
-mam(G 4 Ml)MC 2 )], (13) 
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3.2.4 W boson— neutral Higgs— bottom (top) loop 



There are two contributions, one with a bottom and one with a top quark in the loop (with Hf 
and W interchanged), see Fig. lc. We use the £ = 1 gauge. The WHib loop gives: 



ReSY b CP (WH l b) = - 



V2> 



32tt 2 



lm(X R ) [(3ml ~ 2™|J Im(C ) + m 2 Im(d) + 2m 2 b Im(C 2 ) 
+ Im(5 (m^ + , m 2 w ,m 2 Hi )) - 2lm(B (m 2 ,m 2 b ,m 2 v ))] 
+ m 2 lm(Xt) Im(2C + C 2 )| , (16) 



(17) 



Re5>f P (W^6) = - m t m b i lm (X*) M2Ci + C 2 ) + Im(X fe L ) Im^ - C )] 



where = g HlH+w -S b { R , X b L = g HiH+w _s h { L , and C x = Cx(m 2 t ,m 2 H+ ,m 2 b ,m 2 b ,m 2 w ,m 2 Hi ). 
Analogously, the H\Wt loop gives 



ReSY^^Wt) 



a/2 < 



32tt 2 



Re<5r t CP (IWt) 



32tt 2 



m, m b [lm(Af ) Im^d + C 2 ) + Im(X p ) Im(C7i - C )} , (18) 



Im(X, L ) [(3m 2 - 2m 2 Hi ) Im(C ) + Im(Gi) + 2m t 2 Im(C 2 ) 



+ lm(BQ(m 2 H+ ,m 2 w ,m 2 Hi )) - 2lm(B Q (ml,m 2 ,m 2 v ))] 

+ m\ lm(X R ) Im(2C + C 2 ) } , (19) 



with = g HlH+w - sl' R , X^ = g HlH+w - a]*, and C x = C x (rn%, m 2 H+J m\, m 2 , m\ v m%y) 



t,L 



3.2.5 Ghost— neutral Higgs— bottom (top) loop 

Since the above graphs with a W boson are calculated in the £ = 1 gauge, we also have to 
include the corresponding graphs with W ± — > G ± . These lead to 



Re 5Y b CP (GH l b) = - — m b h b cos f3 Im (Xf ) Im (C ) + m t h t sin /? Im (X 6 R ) Im (d] 

— mbhb cos P Im (X^) Im (C 2 
Re 5Y t CP (GH l b) = — \ m b h t sin (3 Im (X fe L ) Im (C ) + m*/^ cos /5 Im (X 6 L ) Im (Ci 



— m b h t sin /? Im (X^) Im (C 2 



(20) 



(21) 



and 



Re SYf p (HiG t) = - — [ mth cos f3 Im (X?) Im (C ) - m t h b cos /3 Im (X t L ) Im (CV 



Re5Y t CP (H l Gt) 



8vr s 



+ m 6 /i t sin Im (Xf) Im (C7 2 ] 



m t h t sin /3 Im (X t L ) Im (C ) - sin /? Im (Xf ) Im (C x ) 



+ m fe /i fe cos /3 Im (X t L ) Im (C^ 



(23) 



Here, X pL = g HiH+G _ for q = b,t. The C functions are Cx = Cx{fn 2 , m 2 H+ ,ml, m 2 , rn 2 Hi ) 
for q = b and Cx = Cxijn 2 , m 2 H+ ,m 2 , m 2 , m 2 Hi ,m 2 v ) for q = t. 



3.3 Self-energy graphs 
3.3.1 Neutralino— chargino loop 

The self-energy graph with a neutralino and a chargino of Fig. Id gives 



Re6Y b c , P (x° k xt-W) = ± 



9 m h {t) 



){t)y - 8vr 2 V2m 2 H+ m 2 w 



Im ( B (m 2 H+ , m 2 , m 2 

V X k % 



X 



Im (c/j) m -+ (m^+ + m|o - m|+ ) - Im (c u ) m^o (m 2 H+ - m 2 ~o + m 2 + ) (24) 



'Xk 



with c u = F p k O^ + FfcOfc , and cjj = Ff k O^ + i^Og. The overall plus sign is for 5Y b CP , and 
the overall minus sign for 8Y t . 



3.3.2 Stop— sbottom loop 

The graph of Fig. le leads to 
Re5Y b CP (i t ~b J -W) = t 



Sg 2 m h{ t) 



167T 2 m\ 



H+ m w 



/ 9 2 
fflr — 777.7 



Im (GwR&Rfi) Im f B (m^ + , 777^, 777 1) 



(25) 



3.3.3 W± -if, and G ± -ii° loops 

The self-energy graph with jy + and iff is shown in Fig. If. Since we use £ = 1 gauge for the 
in the loop, we have to add the corresponding graph with a ghost, i.e. W ± — > G ± in the 



loop. (The second W propagator can be calculated in the unitary gauge. Hence, no ghost is 
necessary in this case.) The two contributions together give: 

Re^ ( WH° - W) = =F ^ ^L. (2m , _ 2mk _ 3mU x 

o^tt ^2 mi 



H+ m w 



3i (cosp lt + sin p0 2i ) lm(B (m 2 H+ ,m 2 Hi ,m 2 w )) . (26) 



4 Numerical results 

Let us now turn to the numerical analysis. In order not to vary too many parameters, we fix 
part of the parameter space at the electroweak scale by the choice [] 

M 2 = 200 GeV, p = -350 GeV, M d : Mq : M b = 0.85 : 1 : 1.05, 

A t = A b = -500 GeV. (27) 

Moreover, we assume GUT relations for the gaugino mass parameters Mi, M2, M3. In this 
case, the phases of the gaugino sector can be rotated away. Since M , the phase of //, is highly 
constrained by the EDMs of electron and neutron, we take 0^ = 0. The phases relevant to our 
study are thus <pt and the phases of A t and Ab. 

The choice in Eq. Q27|) together with = 490 GeV and tan/5 = 10 gives a mass spectrum 
quite similar to the Snowmass point SPSla [TJJ. Figure shows 5 C7P for this CBjSG clS cL function 



of rriH+ in the range mjy+ = 200- 1400 GeV and various values of <fif The sparticle masses are 

given explicitly in Table 1. 

For mjj+ < Tn^ + m^, 5 CP is very small, O(10~ 3 ) or smaller. The contributions to 5 CP 

come from the diagrams of Figs, la, lc, and If; the diagram of Fig. |T]d only contributes if there 

is a non-zero phase in the chargino/neutralino sector. In the region = 200 — 800 GeV, one 

can distinguish the thresholds of X1X1 a ^ m H+ — 290 GeV, X1X2 a ^ m H+ — 470 GeV, X3X1 at 
5 In the original version of the paper ||l3f there were two mistakes in the analytic expressions [fl4|| , which 
also affected the numerical analysis: for tan/3 = 10 (40) the results for 8 CP in |ll| are reduced to about 30% 
(50%). However, one can easily find scenarios where 5 CP goes up to 10-15%. In particular, for the parameters 
of Eq. (J27|), the resulting plots are very similar to the original ones. 
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490 10 
40 


384 (377) 568 (573) 
379 (377) 571 (573) 


486 522 
435 566 


350 10 
40 


226 (213) 465 (471) 
216 (212) 470 (472) 


340 382 
257 443 



Table 1: Sparticle masses (in GeV) for parameter sets used in the numerical analysis for 
<k = (f )• 

m H + ~ 550 GeV, X2X2 an d xlxt a t — 560 GeV, and X34X2 a ^ m H+ — 730-740 GeV. 

Below the XiXi threshold, 5 CP originates only from the graphs with W and a neutral Higgs 
boson of Figs, lc and If. Here note that the graph with WHf b of Fig. lc always contributes, 
since m t > raw + rrib- 

However, once the H + — > tb channel is open, 5 CP can go up to several per cent. The 
thresholds of H + — > t\b\ at m#+ ~ 860 GeV, and of H + — > t 2 ^2 a t m H+ ~ 1100 GeV are 
clearly visible in Fig. ||. For m H + = 1 TeV, we obtain 6 ~ —3%, —6%, and —8% for (p t = |, 
j, and |, respectively. The dominant contribution comes from the stop-sbottom-gluino loop of 
Fig. |l]b. Also the stop-sbottom-neutralino loop of Fig. lb and the stop-sbottom self-energy of 
Fig. le can give a relevant contribution and should thus be taken into account. The contribution 
of the graphs with x^X or H°W (Fig. |l]a,c,f) exchange can, however, be neglected in this case. 

The relative importance of the various contributions is illustrated in Fig. |||, where we plot the 
form factors R,e5Yf !P and R,e5Y t CP as functions of m#+, for <p t = tt/2 and the other parameters 
as in Fig. To calculate the contributions with neutral Higgs bosons, we have used P, |T!| . 
This is sufficient for our purpose, since we are mainly interested in large CP-violating effects 
that occur for m#+ > m~ tl + mg because of <pt,b- However, once precision measurements of H ± 
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Figure 2: Absolute value of 5 as a function of rriH+, for Mq = 490 GeV and tan/5 = 10. 
The solid, dashed, and dotted lines are for <p t — |, |, and |, respectively, 4>b — 0, the other 
parameters are fixed by Eq. (|2T|). 
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Figure 3: Absolute values of ReSY^ p and ReSY t CP as functions of m H +, for <f) t = n/2 and the 
other parameters as in Fig. |[ The blue lines show the contribution from the x^X exchanges 
of Fig. la; the red lines are those from the tb exchanges of Figs, lb and le; the green lines 
are those from the diagrams with H° and W of Figs, lc and If. The red dashed lines show 
Re 5Yf p due to the tbg loop only. 
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decays become feasible, a more complete calculation of the Hf masses and couplings |15 might 
be used. 

We next lower the stop/sbottom mass scale to Mq = 350 GeV. The resulting masses are 
given in Table 1. Figures (|a and f|b show 5 CP for this case as functions of and tan/3, 
respectively. The threshold behaviour of Fig. [|a is very similar to that of Fig. 0a. The threshold 
of H + — > t\bi is shifted to m#+ ~ 550 GeV, and 5 CP reaches larger values for lighter squarks. 
Even for a small phase <frt, 5 CP can be a few per cent. Figure |4]b shows the tan/3 dependence 
of S CP for m,H+ = 700 GeV and the cases <pt = tt/2, </>& = (full line) and <pt = <frb = k/2 
(dashed line). For completeness, we also show as a dotted line the case of fi = +350 GeV (with 
<pt = 7r/2, <pb = 0). It turns out that the asymmetry has a maximum around tan/3 ~ 10 and 
decreases for larger values of tan/3. In particular, we have 5 CP ~ —12% (—3.5%) for tan/3 = 10 
(40), 4>t = 7r /2, 4>b = 0, and /i = —350 GeV. An additional phase of At, can enhance or reduce 
the asymmetry. For /i < 0, however, its effects in the triangle and self-energy graphs of Fig. lb 
and le compensate each other so that the overall dependence on 0& is small. The dependence 
is larger for \i > 0. 

Here we also note that the branching ratio of H + — > tb increases with tan/3. In the case 
of vanishing phases we have BR(H + tb) ~ 17% (15%) at tan/3 = 10 and 85% (75%) at 
tan/3 = 40 for Mq = 350 (490) GeV and m H+ = 700 (1000) GeV. 

The dependence on <p t is shown explicitly in Fig. |5], where we plot 5 CP as a function of <p t for 
Mq = 350 GeV, m H + = 700 GeV, tan/3 = 10 and 40 and various choices of </V As expected, 
6 shows a ~ sin <j) t dependence. 

Last but not least we relax the GUT relations between the gaugino masses and take rrig 
as a free parameter (keeping, however, the relation between Mi and M 2 ). Figure |B|a shows 
the dependence of 5 CP on the gluino mass, for Mq = 350 GeV, m H + = 700 GeV, 4>t = 7r /2, 
and tan/3 = 10 and 40 (keeping M3 = m~ g real). It is interesting that there is still an effect 
for a large gluino mass: for tan/3 = 10, 5 CP is reduced from about —12% to about —9% for 
600 — * 1200 GeV. Also in the large tan/3 case, 5 CP is decreased by 30% when the gluino 



m 



u 



mass is doubled: from —3.4% to —2.4% for m„ = 600 — > 1200 GeV. A non-zero phase of M3 
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Figure 4: 5 CP for Mq = 350 GeV; in (a) as a function of for tan/3 = 10, <pb = 0, and in 

(b) as a function of tan/3, for m#+ = 700 GeV and (j) t = |. The other parameters are fixed by 
Eq. (0). 
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ft, 




Figure 5: 5 CP function of <j) t , for Mq = 350 GeV, m H + = 700 GeV, tan/3 = 10 in (a) and 
tan/3 = 40 in (b); full lines: 4>b = 0, dashed lines: = dotted lines: /x = 350 GeV, fa = 0. 
The other parameters are fixed by Eq. (P7|). 
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may also have a large effect. This is shown in Fig. ^d, where we plot 5 CP as a function of 03 
for rrtg = \M 3 \ = 565 GeV and t = and w/2. In both cases, the asymmetry can be 0(10%). 
For jj, > 0, the curves are shifted but the order of magnitude does not change. 




Figure 6: a) 5 CP as a function of trig, for Mq = 350 GeV, m#+ = 700 GeV, <p t = 
and tan/5 = 10 and 40. b) 5 CP as a function of 03, the phase of M 3 , for = 565 GeV, 
Mq = 350 GeV, m H+ = 700 GeV, tan/3 = 10, fa = (full line) and & = tt/2 (dashed line). 
06 = 0, the other parameters are fixed by Eq. (P7|). 

5 Conclusions 

We have calculated the difference between the partial rates T (H + — > tb) and V (H~ — > tfe) due 
to CP- violating phases in the MSSM. The resulting decay rate asymmetry 5 , Eq. (1), could 
be measured in a counting experiment. If < + , 5 CP is typically of the order of 
10~ 3 . However, for > + mg , 5 CP can go up to 10- 15%, depending on the phases of 
At and A&, and on tan/3. Also a phase of M3 can have a large effect. 

At the Tevatron, no sensitivity for detecting H ± is expected for a mass > 200 GeV. 
The LHC, on the other hand, has a discovery reach up to ~ 1 TeV, especially if QCD and 
SUSY effects conspire to enhance the cross section. With a luminosity of C = 100 fb _1 , about 
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217 signal events can be expected for pp — > H + ib with S/ \JB = 6.3 (B being the background), 
for m H + ~ 700 GeV and tan/3 = 50 In e + e~ collisions, the dominant production mode is 
e + e~ — > H + H~ . Therefore, one would need a centre-of-mass energy > 2m#+. This would 



certainly be realized at a multi-TeV linear collider such as CLIC [[L7|]. Hence, a CP- violating 
asymmetry 5 CP of a few per cent should be measurable at the LHC or CLIC. 
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A Masses and mixing matrices 



The neutralino mass matrix in the basis of 



(28) 



is: 



/ 



M 



N 



-mz sin 9w cos (3 mz sin 9w sin (3 



Mi 

M 2 m z cos 9w cos (3 —m z cos 9w sin f3 

m z sin 9 W cos (3 m z cos9wcos[3 — p 

^ mz sin9w sin (3 — mz cos 9w sin [3 —ji 



(29) 



with tan (3 = v^/vi. This matrix is diagonalized by the unitary mixing matrix N: 
^Ad^N^ = diag(m^o, m^o, m^o, m^o) , 



(30) 



where m^, k = 1, ...,4, are the (non-negative) masses of the physical neutralino states. 
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The chargino mass matrix is: 

/ M 2 V2m w sm(3 
Mc=\ (31) 
\ V 2 m w cos j3 // 

It is diagonalized by the two unitary matrices U and V: 

WM C V ] = diag(m~+, m~+) , (32) 

where m y i + 2 are the masses of the physical chargino states. 
The mass matrix of the stops in the basis (t L , t R ) is 

M| + m| cos2/3(± - | sin 2 6 W ) + m 2 (A* - //cot /3) m t 



(A - //* cot /3) m t M? + §m| cos 2/5 sin 2 6 W + m 2 



(33) 

Jv[ 2 - is diagonalized by the rotation matrix R l such that R tJf M/~ R l = diag(m 2 i5 m^) and 



1 . We have: 
t R J \ti) 

e 5 cos 6>i — e 5 sin Qz 

r* = i r r , , i (34) 

e 2^tsin0£ e 2^t COS g~ 
Analogously, the mass matrix of the sbottoms in the basis (&l, 6r), 

M| - m| cos 2/3(| - | sin 2 + m 2 . (A* b - /t tan /3) m fe 





R L2 ' 











l 6 



(A - //* tan /3) m fe M~ - lm| cos 2/3 sin 2 W + m l 



(35) 



is diagonalized by the rotation matrix R b such that R bJ< Ai~ R b = diagfm- , m~ ). 

° J b ov 6i' hi' 

In the neutral Higgs sector, we have two CP-even states fa = v / 2Re(7i*) — Vi, i = 1,2, 
and and one CP-odd state a = V2 (- smplm(H\) + cos/3 Im(H 2 )), where Hi and H 2 are the 
two Higgs doublets in the interaction basis. In the basis (</>i, </>2, a), the neutral Higgs mass 
matrix M. 2 H can be written as the well-known tree-level part, which has a block form in this 
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M 



(36) 



basis, plus a general 3x3 matrix containing the loop corrections: 
s 2 P m\ + c 2 P m| — s/3 c(3 (m 2 A + m|) 
-s/3 c/3 (m A + m|) c 2 p" + s 2 /? m| 


where s/3 = sin j3, c(3 = cos j3, etc. In the case of complex parameters, the loop contributions of 
(Ai l H op ) 2 lead to a mixing of the CP-even and CP-odd states. The mass eigenstates then are 



V 



[M l H op \ 



ra\ J 



^3, 



T . 



V a J 



(37) 



The real 3x3 rotation matrix O diagonalizes the mass matrix Ai 2 H , 



T M 2 H O = diag (^n 2 H o,m 2 H o, m 2 H o 



(3* 



with m H o < m H o < m H o. The transformations of the Higgs fields from the interaction basis to 
the mass eigenstate basis are given by 



n\ 

Hi 
Hi 



vi + -= [(O y + z sin P0 3j ) H° - % cos/3G°] , 
v2 

— cos P G~ + sin j3 H~ , 



sinpG + + cospH + , 

v 2 + -^=[(0 2j + i cos /30 3j ) H° 



i sm 



in/5G°] 



(39) 



with the implicit sum over j = 1,2,3, 
V2 



■ R ^ ■ R 

vsmp = m w sm p . 
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v i = v cos p = — m w cos p , v 2 = v sm p — — raw sm P ■ (40) 
For the numerical evaluation of the physical Higgs masses and the rotation matrix O in the 



1-loop effective potential approach [H], we use the program chp.f p"2 



B Interaction Lagrangian 

In this section we give the parts of the interaction Lagrangian that we need for our calculation. 
We start with the interaction of Higgs bosons with quarks and squarks: 

C Hqq = H + i(y b P R + y t P L ) b + H'b (y t P R + y b P L ) t + H? q (sf R P R + sf L P L ) q , (41) 
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C H qq — {G^)^ H + 1* bj + (Gl) i:j H V-ti, 
with i,j = 1, 2, I = 1, 2, 3 and 

^ = 1(1-75), P fi =K 1 + 7 5 )- 
For the H couplings to quarks we have 



y t = h t cos (3, y b = h b sin (3 . 



with 



gm t 



\/2m w sin (3 



hh 



gm b 



a/2 mvi/ cos (3 



where m q is the DR running quark mass. For the H® couplings to quarks we have 



g m q i s j_ • p \ 

\9H iq q + l 9H iq q) ) 

{9H m - % 9u iqq ) > 



2™^ 

£™9 / 5 



(42) 

(43) 

(44) 

(45) 

(46) 
(47) 



2mw 

with gfj iqq and gH iqq given by Eqs. (4.11)-(4.14) in ||. The couplings to squarks are given 
by the matrix 



(48) 



with 



^ h b m b sin (3 + h t m t cos (3 — y/2 g m w sin (3 cos /3 h b (A* b sin (3 + /i cos (3) ^ 
^ /i t (A t cos (3 + /i* sin /3) /itmfc cos /9 + h b m t sin /3 y 



• (49) 



The interactions of charginos and neutralinos are described by 

C Hn = H + ^(FiP R + F^ k P L )xl + H^UFi*P L + F l L k *P R )xt 



C qq ^ 



t(l b yPR + *&P L ) Xfb i + b ^P R + k%P L ) xfK 



+ xt (A;p L + k\-p R )ti: + xf {i%p L + kt;p R )bt* 

9 WPr + b l P L) xUi + xlK k P L + b%P R ) q q* , 

17 



(50) 

(51) 
(52) 



with i,j — 1,2 and k — 1, ...,4. The couplings of H to charginos and neutralinos are 



F R 
r ik 



ik 



-9 

-g 



V iX N u + ^= (iV fe2 + iV fel tan V i2 cos /? , 
- 71 Mi + tan w)U; 2 ] sin £ . 



(53) 
(54) 



The chargino-squark-quark couplings are 



% = -gV jl Rl i + h t V j2 Ei i 
k\j = hb Uj 2 Rii , 



lij — ~9 UjiR^ + h b U j2 R 2i , 
4 = h t V* 2 Rl . 



(55) 
(56) 



The neutralino-squark-quark couplings are 



'■■ik 



9flk R ii + h 9 L k R 2l 



2i ■ 



h Lk R li + 9fRk R 2i > 



(57) 
(58) 



with 



-j^(N k2 + | tan 6 w N kl ) 
^fl tan6 w N* kl , 



h\ k = -ht N* k4 , 



flk 

fnk 

n Lk 



75 (N k2 - \ tan 9 w N kl/ 



V2 

tan^A^fei, 



-hN* k3 . 



(59) 
(60) 
(61) 



Finally, the squark-quark-gluino interaction is given by 



C 



my 



V2g s T« [r(Rl e-&P L - R\, e^P R ) q s q* >t 

+ q s (Ri*e^P R -Rl*e-^P L )~g a qi,t 



(62) 



and the quark interaction with W bosons is 



c„ w = -^ 2 iWft^p L b + w;l^p L t) 



(63) 



We next turn to the interaction of Higgs bosons with W bosons and ghosts. The Lagrangian 
of two Higgs particles and one W boson is given by 



C 



HHW 



9 

V2 



w+ (n\* B i n 2 1 + Hi* &* nfj + w~ (n\ nj* + n 1 2 B i n 



(64) 
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where 

Ad^B = A(d^B) - (d„A)B. (65) 
Using the transformations Eq. fl39|) we get 



r - ■ 9 

i-HHW — 1 ^ 



(g Hj H-w+H° <? H~ + g H]G ~w + H° S t G~ + i G° G~) + h. c. 



(66) 



where 



9H :j H-w+ = — sin /3 Oij + cos j3 2 j + i 3j and (67) 
9h 3 g-w+ = cos/30ij + sin (J 2j ■ (68) 

Moreover, gUjH+w- — 9*HjH-w+ anc ^ QhjG+w- — gH J G~w+- Note that there is no G° W ± H T 
coupling. We further need the couplings of two charged and one neutral Higgs bosons. They 
are derived from the D-term interaction Lagrangian of the Higgs sector, 

C = -\ (D'D' + D 1 D 1 + D 2 D 2 + D 3 D 3 ), where D' and D* are the U(l) Y and ST/(2)j D-terms, 
respectively. In terms of Higgs fields in the interaction basis we have 



P ^ ( ~2 i /'l_/l*'7_/l i n_/2*nj2 n_/l*n_/l nj2*n_/2\ 

J-HHH = — - [9 + g J [rii rti + rii rl\ — rt 2 n 2 — ri 2 n 2 ) 



- £ (n\*nl + nfn 2 2 ) (n\nl* + n\n 2 2 *) . (69) 

The Lagrangian in the mass eigenstate basis is again obtained by applying the transformations 
of Eqs. (|39|) and (|40|). We are only interested in the combinations of (H®, G°) x (H ± , G r± ) x 
(IP, G T ). The couplings to G°, e.g. G°H + H~, G°H + G~, G°G + H~, and G°G + G~, are zero. 
The couplings to Hf are: 



£hhh = gHjH+H- H®H + H + gHjH+G- H®H + G 

+g* HjH+G - H°G+H- + g HjG+G - H°G + G~ , (70) 



with 



g Hj H +H - = 9 -^y~ {[(1 + t 2 w ) c2[3 - 2] c/30 y - [(1 + t 2 w ) c2/5 + 2] s/30 2i } , (71) 

9H jH+G - = 9 -^y~ [rtPWOv + c(30 2j )-i 3j +t 2 w s2/3(c/30 y -s/3 2j )] , (72) 

9 HjG+ G- = - g ^(l+t 2 w )c2(3(c(30 lj -s/30 2j ). (73) 
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Note that only gUjH+G- is complex. In Eqs. (|7"l|) - (ff3|) , we have used the abbreviations s/3 = 

2 

W 



sin/?, s2/? = sin 2/3, c/3 = cos/3, c2/3 = cos 2/3, and tL = tan 2 #w 



C Passarino— Veltman integrals 

Here we give the definition of the Passarino- Veltman one-, two-, and three-point functions || 
in the convention of [ID|. For the general denominators we use the notation 

V° = q 2 - ml and V j = (q + Pj ) 2 - m) . (74) 

Then the loop integrals in D = 4 — e dimensions are as follows: 

Mml) = ^i/A^, (75) 
B (plmlml) = JLjdPq^L;, (76) 
BM,m 2 ,ml) = ^ J d D q^^=pi fl B 1 (p 2 1 ,m 2 Q ,m 2 1 ), (77) 



and 



c » = h\ dDq = Vi » Ci + p ^° 2 ' (79) 

C - — I ' d D a %<lv 

= Q^Coo + pitfPivCn + (pi^p 2v + P2ixPiu)C 12 + p 2 ^P2uC 2 2 , (80) 

where the C"s have (p 2 , (pi — £>2) 2 ,P2> m o> m \i m i) as their arguments. The function B\ can be 
expressed as a combination of the functions A and 5o: 

2p 2 1 B 1 {p 2 1 , m 2 , m 2 ) = A (m 2 ) - A Q {m\) - (p\ -m\ + m 2 ) B (pl, m 2 , m 2 ) . (81) 
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